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Abstract 

We analyze the next-to-leading order SUSY-QCD corrections to the production of 
Higgs particles at hadron colliders in super symmetric extensions of the Standard 
Model. Besides the standard QCD corrections due to gluon exchange and emission, 
genuine supersymmetric corrections due to the virtual exchange of squarks and 
gluinos are present. At both the Tevatron and the LHC, these corrections are found 
to be small in the Higgs-strahlung, Drell-Yan-like Higgs pair production and vector 
boson fusion processes. 



* Heisenberg- Fellow 

Supported by the EU FF Programme under contract FMRX-CT98-0194 (DG 12 - MIHT) 



1 



1. Introduction 



A firm prediction of supersymmetric extensions of the Standard Model (SM) QXJ] is the 
existence of a light scalar Higgs boson. In the Minimal Supersymmetric Standard Model 
(MSSM) the Higgs sector contains a quintet of scalar particles [two CP-even h and H, a 
pseudoscalar A and two charged H ± particles] ||, the Higgs boson h of which should be 
light, with a mass M h < 135 GeV §. If this particle is not found at LEP2 |J, it will be 
produced at the upgraded Tevatron (where a large luminosity, / C ~ 20 fb _1 , is expected) 
§ or at the LHC g |7[, if the MSSM is indeed realized in Nature. 

At hadron colliders, the dominant production mechanism for neutral Higgs particles 
is the (heavy quark) loop induced gluon fusion process, gg —>■ $ with $ = /i, H or A M. 
Since the Higgs particles in the mass range of interest, M$ < 135 GeV, dominantly decay 
into bottom quark pairs, this process is rather difficult to exploit at the Tevatron because 
of the huge QCD background ||. In contrast, at the LHC rare decays of the lightest h 
boson to two photons or decays of the heavy H, A boson to r and \i lepton pairs make 
this process very useful |J. 

Additional Higgs production mechanisms at hadron colliders are provided by: 

a) Higgs-strahlung off W or Z bosons for the CP-even Higgs particles [due to CP- 
invariance, the pseudoscalar A particle does not couple to the massive gauge bosons at 
tree level]: qq -> V* -> $V with $ = h, H and V = W, Z f| . At the Tevatron, the 
process qq' — > hW [with the h boson decaying into bb pairs] develops a cross section of 
the order of a fraction of a picobarn for a SM-like h boson with a mass below ~ 135 GeV, 
making it the most relevant mechanism to study ||. At the LHC, both the bb and 77 
decay modes of the h boson may be exploited ||. 

b) If the heavier H, A, H ± bosons are not too massive, the pair production of two Higgs 



particles in the Drell-Yan type process, qq — > $i<l ) 2 [10-12], might lead to a variety of final 
states [hA, HA, ^h, E ± E, E ± A, H + H-] with reasonable cross sections [in particular 
for Ma ~ Mh ~ M^± < 250 GeV and small values of tan/3, the ratio of the vacuum 
expectation values of the two Higgs doublets] especially at the LHC. Moreover, neutral 
and charged Higgs boson pairs will be produced in gluon fusion gg — > $i$2 |PHl3| . 

c) The production of CP-even Higgs bosons via vector boson fusion, qq — > qqV*V* — > 
||14|| : In the case of a SM-like h boson, this process has a sizeable cross section at 
the LHC. While decays of the Higgs boson into heavy quark pairs are problematic to 
be detected in the jetty environment of the LHC, decays into r lepton pairs make this 



process useful at the LHC [fT5 |. 



In addition to these types of processes, neutral Higgs boson radiation off heavy bottom 
and top quarks [qq, gg — > 66$, tt$] might play an important role in SUSY theories ||16|| . In 



particular, because the couplings of the Higgs boson to b quarks can be strongly enhanced 
for large values of tan/3, Higgs boson production in association with bb pairs can give rise 
to large production rates. 
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It is well known that for processes involving strongly interacting particles, as is the 
case for the ones discussed above, the lowest order cross sections are affected by large 
uncertainties arising from higher order corrections. If the next-to-leading QCD corrections 
to these processes are included, the total cross sections can be defined properly and in a 
reliable way in most of the cases. 

For the standard QCD corrections, the next-to-leading corrections are available for 
most of the Higgs boson production processes. The K-factors [defined as the ratios of the 
next-to-leading order cross sections to the lowest order ones] for Higgs boson production 



via the gluon fusion processes have been calculated a few years ago jTTj; the [two- loop] 
QCD corrections have been found to be significant since they increase the cross sections 
by up to a factor of two. The K-factors for Higgs production in association with a gauge 
boson (a) and for Drell-Yan-like Higgs pair production (b), can be inferred from the Drell- 
Yan production of weak vector bosons and increase the cross section by approximately 



30% ]T8| |. The QCD corrections to gg — > $i$2 are only known in the limit of light Higgs 
bosons compared with the loop-quark mass; they enhance the cross sections by up to a 



factor of two []T3[| . For Higgs boson production in the weak boson fusion process (c), the 
QCD corrections can be derived in the structure function approach from deep-inelastic 
scattering; they turn out to be rather small, enhancing the cross section by about 10% 
[ ]n| . Finally, the full QCD corrections to the associated Higgs production with heavy 
quarks are not yet available; they are only known in the limit of light Higgs particles 
compared with the heavy quark mass [M which is only applicable to tih production; in 
this limit the QCD corrections increase the cross section by about 20-60%. 

Besides these standard QCD corrections, additional SUSY-QCD corrections must be 
taken into account in SUSY theories; the SUSY partners of quarks and gluons, the squarks 
and gluinos, can be exchanged in the loops and contribute to the next-to-leading order 
total cross sections. In the case of the gluon fusion process, the QCD corrections to the 
squark loop contributions have been calculated in the limit of light Higgs bosons and 
heavy gluinos; the K-factors were found to be of about the same size as the ones for the 
quark loops pl| . The SUSY-QCD corrections to the other production processes are not 
yet available. 

In this paper, we calculate the SUSY-QCD corrections to the Higgs production cross 
sections for Higgs-strahlung, Drell-Yan-like Higgs pair production and weak boson fu- 
sion processes. These corrections originate from qqV one-loop vertex corrections, where 
squarks of the first two generations and gluinos are exchanged, and the corresponding 
quark self-energy counterterms, Fig. |I[ These corrections turn out to be small for rea- 
sonably large values of the squark and gluino masses, altering the cross section by a few 
percent only. For heavier masses, squarks and gluinos decouple and the SUSY-QCD cor- 
rections become tiny. 

Note that genuine SUSY-QCD corrections are also present in the decays of the Higgs 
bosons into heavy quark pairs fl22|. They are much larger than for the production in the 
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Figure 1: Generic diagrams contributing to the SUSY-QCD corrections to the qqV vertex 
[V = i,Z,W] at NLO. 



case of large mixing in the sbottom sector, being of the order of a few ten per cent, and 
decouple only slowly for large squark and gluino masses. They are being included in the 
program HDECAY |23j which calculates the decay widths and branching ratios of the 
MSSM Higgs particles. 

2. QCD corrections to the production processes 

2.1 Higgs-strahlung and Drell— Yan-like pair production processes 

At hadron colliders, the lowest order partonic cross section for the Higgs-strahlung pro- 
cesses, qq — > V$ with V = W, Z and $ = h,H, is given by J?| [| 

A , G 2 F My 2 . 2 2^1/2/^2 ~MM&,M£;s) + VlMlls 

*to(w - VS) = -^jdlvvH + O^iMl Ml- s) K ^J^TJ^ ~ 0) 

with the couplings g$vv — sin(/3 — a) or cos(/9 — a) for h and if respectively; s is 
the partonic cm. energy and A the usual two-body phase space function X(x, y; z) = 
(1 — x/z — y/z) 2 — Axy/z 2 . v q , a q are the vector and axial- vector couplings of the quark q 
to vector bosons and are given by v q = 2I q — 4e g s^, a q = 2I| for V = Z [e q is the electric 
charge, I q the weak isospin of the quark and sly = 1 — c\y = sin 2 6 W ] and v q = a q = y/2 
for V = W. 

The partonic cross section for Drell- Yan-like Higgs pair production qq —> $i$2) where 
at least one of the Higgs bosons is neutral, is given by [[LT], [TTJ : 

where the couplings are given by: g ZhA = g WhH + = cos(/3 - a), g ZHA = g WH H+ = 
sin(/3 — a) and gwH+A = 1- For charged Higgs boson pair production there is an s- 
channel photon exchange in addition to the Z-channel diagram, and the cross section is 
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given by (TO, |T2 



Tra 2 (s) 
9s 



2e q v„v H + 



+ 



(v 2 + a 2 )v 2 



q) v H+ 



16s 2 w c 2 w (l - Mljs) 64s 4 w c 4 w (l - M 2 /sf 



(3) 



with v H + = — 2 + 4sfy. 

The hadronic cross sections can be obtained from convoluting eqs. (1-3) with the 
corresponding (anti) quark densities of the protons 



TS 



(4) 



where r = (My + M$) 2 /s for the Higgs-strahlung and r = (M $1 + M$ 2 ) 2 /s for the pair 
production processes, with s being the total hadronic cm. energy squared. 

The standard QCD corrections, with virtual gluon exchange, gluon emission and quark 
emission, are identical to the corresponding corrections to the Drell-Yan process flT8| . p"-f . 
They modify the lowest order cross section in the following way |7|, [18 
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with the coefficient functions 124 
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uj qg (z) = 
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where M denotes the factorization scale, /i the renormalization scale and P qq , P qg the 
well-known Altarelli-Parisi splitting functions, which are given by j2S] 



PqqW 
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l-z) 



P qg {z) = \{z 2 + {l-z) 2 } 



(7) 



The index + denotes the usual distribution F + (z) = F(z) — 5(1 — z) f dz'F(z'). 

Including the vertex correction due to the squark-gluino exchange diagram and the 
corresponding self-energy counterterm, the lowest order partonic cross section in eq. (4) 
will be shifted by 

2 a s (jj) 



olo — * 0"LO 



1 + 



3 7T 



-3?eC(s, m q , m g ) 



(8) 
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For degenerate unmixed squarks [as is approximately the case for the first two generation 
squarks], the expression of the factor C is simply given by 



»1 rl 

C(s,mq,mg) = 2 I xdx / dylog- 



mj + (m| — m~)x 



(9) 



sx 2 y{l — y) + (m| — m~)x + m~ — ie 

the expression of 



In terms of the Passarino-Veltman scalar functions Aq, B and Cq 
the function C reads 

2 

C(s,mq,mg) = 2m?C (s, 0, m,-, m,-, m^) - 7 (m^- - m^)C + (s, 0, m,-, m,-, m^) 



1 + B (s, uig, Trig) - 2Si(0, m g , m q 



(10) 



with 

C + (s, 0,m 1 ,m 1 ,m 2 ) 
B 1 (0,m 1 ,m 2 ) 

where 



B Q (s, mi, mi) - B o (0, m 1 ,m 2 ) + (m 2 , - Tn\)C (s, 0, m l ,m l ,m 2 ) 
i B (0,mi,m 2 ) + (ml — ml)B' (0,Trii,Tri 2 ) (11) 



d 

B' (s,mi,m 2 ) = —B (s,mi,m 2 ) 



:i2) 



It should be noted that these results agree with the corresponding SUSY-QCD corrections 
to slepton pair production at hadron colliders |27[ . 



2.2 The vector boson fusion processes 

The differential leading order partonic cross section for the vector boson fusion process, 
qq — > qqV*V* — > gg$ with V = W, Z and $ = h,H [with the couplings g<s>vv given before] 
can be cast into the form JF], [L4j] 



da 



LO 



1 y/2G 3 F M*q>qigl vv 
4[ ?1 2 - M 2 ] 2 [g 2 2 - MIY 

Fi(xi,M 2 )F 2 (x 2 ,M 2 ) 

P2Q2 

F 2 (xi,M 2 )Fi(x 2 ,M 2 ) 
Pi?i 

F 2 { Xl ,M 2 )F 2 {x 2 ,M 2 ) 

(Pigi)(P 2 g 2 ) 



F 3 (x 1 ,M 2 )F 3 (a; 2 ,M 2 
2(P iqi )(P 2 q 2 ) 



F 1 {x 1 ,M 2 )F 1 (x 2 ,M 2 



2 + 



(glg2) ! 

Sid 



(Ag 2 ) 2 A/f2j _ 1 / p P 2 g 2 N 
— = M P + -J P 2 gi 2-gig 2 



- M P + — \ Piq 2 — qiq 2 

?2 \ Ql j 



P1P2 



(Pigi)(P 2 gi) (P 2 g 2 )(Pig 2 ) 



9? 



4 



(Pigi)(P 2 g 2 )(gig 2 ) 



[(P 1 P 2 )(g 1 g 2 ) - (Pig 2 )(P 2 gi)] G^ob, 



dPS, 



(13) 



where dPS 3 denotes the three-particle phase space of the final-state particles, M P the 
proton mass, P\p the proton momenta and the momenta of the virtual vector bosons 
V*. The functions F^x, M 2 ) {i = 1,2,3) are the usual structure functions from deep- 
inelastic scattering processes at the factorization scale M: 

Fx(x,M 2 ) = $>J + a;)[g(x,M 2 ) + g(x,M 2 )] 
<? 

F 2 (x,M 2 ) = 2xJ2(v 2 q +a 2 q )[q(x>M 2 ) + q(x,M 2 )} 
g 

F 3 (x,M 2 ) = 4Y,v q a q [-q(x,M 2 ) + q(x,M 2 )} (14) 
</ 

where v q and a q have been defined previously. 

In the past the standard QCD corrections have been calculated within the structure 
function approach fl9 |. Since at lowest order, the proton remnants are color singlets, at 



NLO no color will be exchanged between the first and the second incoming (outgoing) 
quark line and hence the QCD corrections only consist of the well-known corrections to 
the structure functions Fi(x, M 2 ) (i = 1,2,3). The final result for the QCD-corrected 
cross section can be obtained from the replacements 



Fi(x, M 2 ) -> Fi(x, M 2 ) + AFi(x, M 2 , Q 2 " 



(i = 1,2,3) 
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where z = x/y and the Altarelli-Parisi splitting functions P gg , P gff are as given before. We 
introduced the notation T>iiz) = (log*(l — z)/(l — z)) + (i = 0, 1). The physical scale Q is 



given by Q 



-q~ for x 



X i 



1, 2). These expressions have to be inserted in eq. (|13|) 



and the full result expanded up to NLO. The typical renormalization and factorization 
scales are fixed by the corresponding vector-boson momentum transfer /i 2 = M 2 = —qf 
for x = Xi (i = 1, 2). 

The inclusion of the (two) vertex corrections due to the squark-gluino vertex diagrams 
and the corresponding self-energy counterterms, assuming that the quarks are massless 
and the squarks degenerate in mass, can be performed by shifting the functions F^Xj, M 2 ) 
in eq. (13) by (i = 1, . . . , 3 and j = 1, 2) 



F t (x v M 2 ) 



Fi(x v M 2 ) 



1 + |— jj^-3fceC(gJ, rrig, m- g 



(19) 



where the function C(s,mq,rrig) has been given in eq. fllCf) . 



3. Numerical Results 

We will perform our numerical analysis for the light scalar Higgs boson h in the decoupling 
limit of large pseudoscalar masses, Ma ~ 1 TeV. In this case the light h boson couplings to 
standard particles approach the SM values. In our analysis we will investigate the Higgs 
strahlung process qq — > hV and the vector boson fusion mechanism qq — > qqV*V* — ► qqh 
[V = W,Z). 

In Fig. 2, we show the cross sections for these processes at LO and at NLO with only 
the standard QCD corrections included, as a function of the h boson mass [which can be 
made varied by varying the parameter tan /3] for Tevatron and LHC energies. The NLO 
(LO) cross sections are convoluted with CTEQ4M (CTEQ4L) parton densities and 
NLO (LO) strong couplings a s . As can be inferred from the figure, the standard QCD 
corrections increase the Higgs-strahlung cross sections by about 30% (40%) at the LHC 
(upgraded Tevatron) and the fusion process by about 10% (5%) at the LHC (upgraded 
Tevatron) . 

For the SUSY-QCD corrections, we evaluated the Higgs mass for tan /3 = 30, Ma = 1 
TeV and vanishing mixing in the stop sector; this yields a value Mh = 112.6 GeV for 
the light scalar Higgs mass. For the sake of simplicity we decompose the K factors 
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Figure 2: Total LO and NLO cross sections of Higgs boson production via Higgs-strahlung 
qq — > h + W/Z and vector boson fusion qq — > qqV*V* — > qqh [V = W, Z] at the LHC 
(left) and the Tevatron (right) in the decoupling limit. 
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Figure 3: Relative corrections due to virtual squark and gluino exchange diagrams to 
Higgs boson production via Higgs-strahlung qq — > h + W/Z and vector boson fusion qq — > 
qqV*V* -> qqh [V = W,Z] at the LHC (left) and the Tevatron (right). 
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K = pnlo/clo into the usual QCD part Kq C d and the additional SUSY correction 
Ssusy' K = Kqcd + &SUSY- The additional SUSY-QCD corrections Ssusy are presented 
in Fig. |^ as a function of a common squark mass for a fixed gluino mass rrig = 200 GeV 
[for the sake of simplicity we kept the stop mass fixed for the determination of the Higgs 
mass Mh and varied the loop-squark mass independently]. 

The SUSY-QCD corrections increase the Higgs-strahlung cross sections by less than 
1.5%, while they decrease the vector boson fusion cross section by less than 0.5%. The 
maximal shifts are obtained for small values of the squark masses of about 100 GeV, which 
are already ruled out by present Tevatron analyses |29| ; for more reasonable values of these 



masses, the corrections are even smaller. Thus, the additional SUSY-QCD corrections, 
which are of similar size at the LHC and the Tevatron, turn out to be very tiny. For large 
squark/gluino masses they become even smaller due to the decoupling of these particles, 
as can be inferred from the upper squark mass range in Fig. ^. 



4. Conclusions 

In this letter we have determined the SUSY-QCD corrections to Higgs boson produc- 
tion via Higgs-strahlung, Drell-Yan-like Higgs pair production and vector boson fusion 
at hadron colliders. The additional SUSY corrections originating from the exchange of 
virtual squarks and gluinos range at the per cent level and are thus rather small. This 
analysis completes the theoretical calculation of the NLO production cross sections of 
these processes in the framework of supersymmetric extensions of the Standard Model. 
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